Frustration-Induced Two Dimensional Quantum Disordered Phase in Piperazinium 

Hexachlorodicuprate 
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Piperazinium Hexaclilorodicuprate (PHCC) is siiown to be a frustrated quasi-two-dimensional quan- 
tum Heisenberg antiferromagnet with a gapped spectrum. Zero-field inelastic neutron scattering and 
susceptibility and specific heat measurements as a function of applied magnetic field are presented. 
At r = 1.5 K, the magnetic excitation spectrum is dominated by a single propagating mode with a 
gap, A = 1 meV, and bandwidth of « 1.8 meV in the {hQl) plane. The mode has no dispersion along 
the b* direction indicating that neighboring a-c planes of the triclinic structure are magnetically 
decoupled. The heat capacity shows a reduction of the gap as a function of applied magnetic field in 
agreement with a singlet-triplet excitation spectrum. A field-induced ordered phase is observed in 
heat capacity and magnetic susceptibility measurements for magnetic fields greater than Hci ~ 7.5 
Tesla. Analysis of the neutron scattering data reveals the important exchange interactions and 
indicates that some of these are highly frustrated. 



I. INTRODUCTION 

Among physical systems that display collective macro- 
scopic quantum phenomena, interacting spin systems are 
perhaps the most experimentally informative. Many 
qualitatively different model systems are available, and 
there are numerous experimental tools providing ac- 
cess to spatial and temporal correlations on microscopic, 
mesoscopic, and macroscopic scales. Owing to their rela- 
tive simplicity, much attention has been devoted to one- 
dimeiisional systems, such as the spin S. 1 Haldane 
chain,EJ the even-leg S=l/2 spin laddergj3 □ and dimer- 
ized (alternating bond) S=l/2 chains.ErLl However, it has 
long been a quest for theorists and experimentalists alike 
to find analogous cooperative singlet ground state mag- 
netism in higher dimensions. 

For the two-dimensional S=l/2 quantum Heisenberg 
antiferromagnet (2D QHAFM), a non-trivial quantum- 
disordered spin-gap phase with a singlet ground state in 
the vicinity of a quantum-critical point has been stud- 
ied extensively by theorists in relation to the roagnetic 
properties of layered cuprate superconductorsB'Ll Sev- 
eral models were suggested to fall into this quantum- 
disordeced phase, ranging from simple quantum dimer 
modelaij, and a valence bond crystal favored by frustra- 
tion on the 2D souate latticalil'Ej, to the quantum kagome 
antiferromagnefll3il3 and the long-spught spin-liquid res- 
onating valence bond (RVB) state. E£lEZl While there are 
numerous materials which, like the layered cuprates, fall 
into the renormalized-classical region of the phase dia- 
gram of the 2D QHAFM, realizations of quantum disor- 
dered 2D spin systems close to the quantum critical point 
are scarce. 

There are higher dimensional coupled spin dimer 
systems such as the three-dimensional dimer networks 



Cs3Cr2Br9[l-find KCuClsjii and the quasi-2D material 
BaCuSi206H3 However, as with strgngly dimerized linear 
chains such as Cu(N03)2 • 2.5D20ja correlations in these 
systems are dominated by isolated spin pairs. For differ- 
ent reasons, this is also the case in the newly discovered, 
two-dimensional frustrated spin system SrCu2(B03)2.tll 
In that material, triplet excitations are localized, because 
of the frustrating symmetry of the interactions.EHlHere we 
report the discovery of a frustration-induced spin-singlet 
phase with a gapped spectrum in the metalo-organic 
compound piperazinium hexachlorodicuprate (PHCC). 
The important distinction from other D > 1 singlet sys- 
tems is that there are uniform extensive paths of interac- 
tion spanning the quasi-2D plane and the strongly cor- 
related spin cluster in PHCC involves several spin pairs. 
Even so there is a spin gap in the excitation spectrum, 
and no magnetic phase transition in the absence of a 
magnetic field. 

These results were obtained through an extensive char- 
acterization of PHCC using magnetic susceptibility, spe- 
cific heat, and inelastic neutron scattering measurements. 
Our findings correct previous reports thai-spggest that 
PHCC is an alternating spin chain system.o In agree- 
ment with previous thermodynamic data, we find a 
triplet mode with a spin gap A = 1 meV, and a band- 
width of 1.8 mcV. However, rather than the dispersion 
being confined to the putative chain direction, c, we 
observe dispersion throughout the a-c plane. In addi- 
tion, our analysis of the wavevector dependent scatter- 
ing intensity shows that no fewer than five spin pairs 
are strongly correlated. Two of these provide positive 
contributions to the ground state energy, indicating the 
presence of frustration. We also present evidence for 
a field-induced magnetic phase transition that provides 
an interesting example of quantum critical behavior in a 



frustrated quasi-two-dimensional spin system. 



II. PIPERAZINIUM HEXACHLORODICUPRATE 

Piperazinium hexachlorodicuprate (PHCC), 
(C4Hi2N2)Cu2Cl6, lias a triclinic crystal structure with 
space group PI, and room temperature lattice constants 
a = 7.984(4) A, b = 7.054(4) A, c = 6.104(3) A, and 
a = 111.23(8)°, (3 = 99.95(9ji. and 7 = 81.26(7)°. (We 
use the nomenclature of Ref.E3.) The lattice parameters 
at temperature T = 1.5 K were found to be a = 7.82(2), 
& = 6.7(3), and c = 5.8(2) A. 

The crystal structure of PHCC is depicted in Fig. |^. 
Figure |l|(a) is a perspective drawing of the crystal struc- 
ture as viewed along the a axis. Copper-chlorine layers 
in the a-c plane are clearly visible. These are well sep- 
arated from one another by the piperazinium dication 
rings. By rotating the structure in Fig. ^[a) about the c 
axis, one obtains the view shown in Fig. y(b), which de- 
picts a single copper-chlorine plane. The Cu^+ ions have 
a distorted 4-1-1 coordination with their neighboring chlo- 
rine atoms. The four short Cu-Cl bonds, with an aver- 
age bondlength of 2.30 A, are approximately co-planar, 
with the long 2.62 A Cu-Cl bond nearly perpendicular 
to this plane. With this coordination, the Cu^+ spin 
density lies predominantly in the four-bond plane, and 
therefore Cu-Cl-Cu superexchange interactions that in- 
volve only the short Cu-Cl bonds should be stronger than 
those involving the long bonds. On this basis, the mag- 
netic Cu-Cu interaction indicated by Bond J_in Fig. |l|(c) 
was predicted to be stronger-thaji Bond 2,Ej leading to 
the alternating chain model.cJ^Ej However, our neutron 
scattering results show that while the spins interacting 
through Bond 1 are indeed the most strongly correlated, 
Cu-Cu couplings not previously considered, such as those 
due to the halide-halide contacts shown as dotted lines in 
Fig. 0(b), result in dispersion along the a direction that 
is stronger than that along c. As we will show, PHCC 
contains a two-dimensional network of magnetic interac- 
tions in the a-c plane. The Cu-Cu interactions that we 
will consider in our analysis are numbered 1 through 8 
in Fig. |l](c). Our results pertaining to these bonds are 
listed in Table |. 



III. EXPERIMENTAL TECHNIQUES 

Powder samples of PHCC were prepared by rapid cool- 
ing from 50 °C to °C of 38% hydrochloric acid solutions 
containing piperazinium dihjidrochloride and copper(H) 
chloride in a 1:4 molar ratioEj Single crystals were grown 
by slowly reducing the temperature of similar saturated 
solutions from 50 °C to 18 °C over 75 hrs. Seed crys- 
tals obtained in this manner were suspended in satu- 
rated solutions and grown further. The crystals typically 
grow as dark red tablets with (100) faces providing the 



largest facets. Deuterated single crystals for the inelastic 
neutron scattering measurements were produced in this 
manner from commercially available piperazinium(d8) di- 
hydrochloride and anhydrous copper chloride(II), using 
35% DCl in D2O as the solvent. The crystals were 89(1)% 
deuterated, as determined by neutron activation analysis. 

Low field DC magnetic susceptibility measurements 
were performed on a powder sample of mass m = 16 
mg in the temperature range 1.7 K < T < 270 K using a 
SQUID magnetometer. AC susceptibility measurements 
using a balanced-coil susceptometer were performed in 
the temperature range 0.125 K < T < 8 K on single 
crystals of typical mass m = 20 mg, in DC fields up 
to 9 Tesla oriented parallel to the c axis. Specific heat 
measurements were performed on a single crystal of mass 
m = 4 mg for 0.125 K < T <-2.5 K in fields up to 9 Tesla 
using relaxation calorimetryo with the field oriented par- 
allel to the a* axis. 

The sample used for inelastic neutron scattering mea- 
surements consisted of three deuterated single crystals 
with a total mass of m = 3.32 grams coaligned within 
0.8°. Measurements in the (Okl) and {hOl) scattering 
planes were performed on the SPINS cold neutron triple 
axis spectrometer at the National Institute of Standards 
and Technology (NIST) Center for Neutron Research in 
Gaithersburg, MD. The horizontal beam coUimation be- 
fore the sample was 50' /ki (A^^) - 80' for the (Okl) mea- 
surements, and 50' /ki (A^^) - 76' for the (hOl) measure- 
ments. A liquid nitrogen cooled BeO filter was placed 
after the sample, and data were collected at fixed final 
energy Ef — 3.7 meV. A horizontally focusing pyrolytic 
graphite (PG(002)) analyzer with acceptance angles of 
2.7° X 7.2° in the horizontal and vertical plane respec- 
tively was used to increase the count rate at the ex- 
pense of broadening the instrumental wavevector reso- 
lution perpendicular to the scattered neutron wavevec- 
tor, kj. At energy transfer huj = 0, the energy resolu- 
tion was dhuj =. 0.13 meV, and the instrumental reso- 
lution ellipsoidcZl projected on the scattering plane had 
Full Width at Half Maximum (FWHM) principal axes 
SQi = 0.043 A^i and 5Q2 = 0.076 A-\ 

Measurements in the (hkO) plane were performed on 
the BT2 thermal neutron triple axis spectrometer at 
NIST. Horizontal coUimations of 60' - 20' - 20' - 60' 
were used, and data were collected at Ef = 13.7 meV 
with a PC filter before the analyzer. This gave typical 
FWHM resolutions at ?ia; = of SHw = 0.84 meV, and 
SQ\\ = 0.02A"i and SQj_ = 0.09A~^ parallel and perpen- 
dicular to Q, respectively. All data were converted to the 
normalized scattering intenaity /(Q, hto) following a pro- 
cedure detailed elsewhere,c3 using the incoherent elastic 
scattering of the sample as measured in each experimen- 
tal configuration. 



IV. EXPERIMENTAL RESULTS 



A. Magnetic Susceptibility 

The DC magnetic susceptibihty x(T) = M/H for a 
powder sample of PHCC is shown in Fig. ^. The data 
were taken in an applied field H — 50 Oe, and the mag- 
netization M was found to be linear in field up to if = 1 
Tesla at T = 1.8 K. As shown in the inset of Fig. |, x{T) 
rises with decreasing T to a rounded maximum at T « 12 
K, foUowedpbjy a rapid decrease, as has been previously 
observcd.OEj The low temperature behavior of x(T) is 
shown in the main panel of Fig. ^, which also includes 
single- crystal AC susceptibility data measured down to 
T = 0.125 K. Apart from a paramagnetic background at 
the lowest temperatures that is attributable to residual 
impurities, these data show the exponentially activated 
dependence of x{T) characteristic of a gapped Heisenberg 
antiferromagnet . 

The AC magnetic susceptibility for a single crystal as 
a function of applied magnetic field at fixed tempera- 
ture is shown in Fig. |^. The sharp feature in the data 
indicates an abrupt change in the magnetic density of 
states of PHCC as the spin gap closes at the critical field 
Hci w 7.5 T at T = 0.125 K. This feature moves to larger 
magnetic field as the temperature is increased, and even- 
tually becomes smeared out due to thermal population 
of multiple energy levels above the spin gap. 



space that were probed. At every wavevector measured, 
the spectrum is dominated by a single, resolution-limited 
mode. Representative energy scans are shown in Figs. |^, 

0, and H They illustrate dispersion of this mode along 
the h and I directions and the absence of dispersion as 
wavevector transfer varies along the k direction. The 
mode disappears upon heating to T=50 K, as shown in 
Fig. 0(a), demonstrating that the observed scattering is 
magnetic in origin. The mode has its minimum energy 
(spin gap) A = 1 meV at Q = (0.5,0,1.5) [Fig. |(a)]. 
The full dispersion of this mode is summarized in Figs. |^, 
|To| and ^ Each point in these figures was obtained by 
fitting a scan such as those in Figs. ||, ^, and |8|to a Gaus- 
sian lineshape. From these figures, the two-dimensional 
nature of the magnetic interactions of PHCC becomes 
clear. Figs. ||(d) and ^(d) show significant dispersion 
along lines in reciprocal space that intersect the two- 
dimensional magnetic zone center at Q = (0.5,0,1.5), 
with somewhat stronger dispersion along h than along 

1. Along the edges of the magnetic zone, additional dis- 
persion with h is seen in Figs, ^(a)-(c), but only min- 
imal variation in peak position was observed with / in 
Figs. |l](a)-(c). From Fig. an upper bound of 0.2(2) 
meV can be set on variation of the mode energy with k. 
The solid curves in Figs. |6[]ll| result from a fit using a 
two-dimensional model for the dynamic spin correlation 
function as described below. 



B. Specific Heat 

Figure || depicts the heat capacity as a function of tem- 
perature measured for different fixed magnetic fields up 
to iJ = 9 T. The lattice contribution to the heat capac- 
ity has not been subtracted from the data. The expo- 
nentially activated heat capacity, which is apparent from 
the inset to Fig. ^ for fields H < 8 Tesla, is further di- 
rect evidence for a spin gap in PHCC. As H is increased, 
there is an increase in the specific heat at lower temper- 
atures and the activation energy decreases, effects that 
indicate that the spin gap is decreasing. For iJ > 8 T, 
there is a singularity in the temperature dependence of 
the specific heat indicative of a phase transition. This im- 
plies that above Hci PHCC enters a 3D-ordered phase, 
as is typical of quasi-lD and quasi-2D gapped quantum 
antiferromagnets. Both the transition temperature and 
the the spin entropy removed from the system via the 
transition increase with H. 



C. Inelastic Neutron Scattering 

To map out the wavevector dependence of the mag- 
netic excitation spectrum of PHCC, energy transfer scans 
at constant Q were performed at 93 different locations 
throughout the reciprocal space of PHCC in zero mag- 
netic field at T = 1.5 K. Figure || depicts the three recip- 
rocal lattice planes and the particular points in reciprocal 



V. DISCUSSION 

A. Magnetic Susceptibility 

The zero field susceptibility data yield significant infor- 
mation about the microscopic spin Hamiltonian. First, a 
measure of the size of the spin gap may be obtained. In D 
dimensions, the asymptotic low-temperature susceptibiJ-i 
ity of a gapped spin system with quadratic dispersionc3 
is proportional to 

X(r) OC T(^/2)-lg-A/fc,T^ (1) 

The solid line in the main panel of Fig. ^ is a fit to the 
two-dimensional form of Eq. ^ including a diamagnetic 
background and a low-temperature Curie tail. The fit 
was restricted to T < 5 K, and gave a spin gap A = 
1.2(1) meV, in good agreement with that observed in the 
neutron scattering data. 

From the high temperature Curie- Weiss behavior of 
x(r) we obtain a value for the sum of all relevant ex- 
change constants in the system 

Jo = $]Jd, (2) 

d 

where {d} is the set of bonds connecting a spin to its 
neighbors. Fitting the sum of a diamagnetic and a Curie- 
Weiss term to our x(r) data in the temperature range 



50 K < r < 250 K we extracted the Curie- Weiss temper- 
ature, Qcw, the average 5— factor and the diamagnetic 
term resuhing from the sample and sample holder. The 
fit is shown as a solid line in the inset to Fig. ||. The 
Curie- Weiss temperature &cw = —19.0(7) K is related 
to Jo through Qcw = ^^11^.0 Solving for the sum of 
the independent exchange constants, we find Jq = 3.3(1) 
meV while gav — 2.2(3). A similar estimate for Jq is 
obtainedj-by fitting x(^) to an alternating spin chain 
model.E3Cil This reflects the insensitivity of the zero-field 
susceptibility to the geometry of the magnetic interac- 
tions. 



B. Specific Heat 

The change in the activated temperature dependence 
of the specific heat with increasing magnetic field for H < 
8 T gives a measure of the reduction of the spin gap as 
the lower critical field is approached. The asymptotic 
low-temperature form for the specific heat, derived under 
the same conditions as for Eq. |^ is given by 

Cp{T) OC r(^/2)~2g-A/fc«T_ (3) 

The solid lines in Fig. ^(b) are fits to this form for D — 
2, with the addition of a term proportional to /T^ 
to account for nuclear spin contributions at the lowest 
temperatures. The derived gap values are shown versus 
the applied field in Fig. ^(c). The approximately linear 
reduction of the spin gap with field is consistent with the 
Zeeman splitting of a triplet of excited states. 

The positions of the peaks observed in the temperature 
dependent specific heat for iJ > 8 T outline the lower 
edge of the ordered region of the H-T phase diagram. 
This phase boundary is shown in the inset to Fig. |^. 
A rough estimate for the upper critical field is Hc2 — 
Jo/gavlJ^B = 26(3)T. Clearly, further work is required to 
map out the full extent of the ordered phase and also 
to determine the structure and symmetry of the order 
parameter. The locations of peaks in the field-dependent, 
single crystal magnetic susceptibility data are also shown 
in the inset to Fig. 0^ Taking into account the anisotropic 
g-tensor of PHCCo, these data are entirely consistent 
with the specific heat data. 

C. Inelastic Neutron Scattering 

The magnetic contribution, /m(Q,?iaj), to the nor- 
malized neutron scattering intensity /(Q, ?iu;) is related 
to the pduiiamic spin correlation function iS"^(Q,w) as 
foUowsEMll 

iMM) = / o^^g'McJ'7^Q^^(Q-Q',t^-cJ') 



Here F{Q) is the magnetic form factor of the Cu 
ion, andJZQ^ is the normalized instrumental resolution 
functionBJ We assume that the magnetic Hamiltonian 
is isotropic in spin space, and thus since no symmetry- 
breaking ordering transition occurs in PHCC at zero 
field, the dynamic spin correlation function should also 
be isotropic. This is supported by previous susceptibility 
measurements, which indicate that exchange anisotropy 
in PHCC is very smallcfl, and by the absence of any ob- 
servable anisotropy splitting of the excitations observed 
by neutron scattering. In this case the different spin po- 
larizations in Eq. ^ can be summed to give 

E {^-P - Q'o.Q'p) S"^iQ',^') = 25""(Q',a;'). (5) 

Owing to the low symmetry of PHCC, there are a large 
number of potential Cu-Cu magnetic interactions within 
the a — c plane. Information about the relative impor- 
tance of spin pair correlations to the ground state energy 
can be obtained in a mo^l independent way through the 
first moment sum rule.tj For an isotropic spin system, 
this is 

/oo 
-00 

= -^E'^d(So-Sd)(l-cosQ.d), (6) 

d 

where Jd is the exchange strength and (Sq • Sd) is the 
two-spin correlation function for the spin-pair with bond 
vector d. We write the Hamiltonian in the form 

'H= 2^^J<i^r -Sr+d, (7) 
r,d 

where the index {r} runs over all spins. The ground state 
energy per spin Eq is closely related to the first moment, 
being simply 

i?o = ^5I^d(So-Sd). (8) 

d 

The first frequency moment of the dynamic correlation 
function 1i{uj)q measured in PHCC is shown in Fig. ^ 
and Figs. p"2[]l^ . The data points were obtained from the 
same Gaussian fits that were used to determine the mode 
energy. 

Unlike what is observed in strongly dimerized spin gap 
systems, the variation of ?i(w)q with Q in PHCC cannot 
be accounted for by a single spin pair. Instead, several 
crystallographically distinct spin pairs are strongly cor- 
related and contribute significantly to the ground state 
energy. This renders the system considerably more in- 
^tricate to describe theoretically. For example, the RPA 
theory based on a dimerized ground state that accounts 
for strongli2.|-(iime«ized systems in one, two, and three 
dimensionsQullSiHj is not adequate here since a single 



spin pair controls the first moment of S°'"{Q,uj) in that 
model. 

In the absence of an adequate theory of gapped spin 
systems for D > 1 that goes beyond a dimer-based ex- 
pansion, we parametrize the measured dispersion i?(Q) 
in PHCC with the following phcnomcnological expres- 
sion, consistent with Bloch's theorem: 

£'(Q) = {Bo + Bh cos(27r/i) + Bi cos(27rZ) 

+ Bhi[cos{2Tr{h + I)) + cos(27r(/i - /))] 

+ B2h cos(47r/i) + B21 cos(47r0)^/^. (9) 

We note that simpler dispersion relations of the form 
E{Q) — Aq + cos(Q • Ri), which provide an effec- 

tive description of related systems with weaker disper- 
sion such as CuHpCl,E£l do not provide a good descrip- 
tion of the dispersion in PHCC. Because most of the ob- 
served magnetic intensity comes in the form of resolution- 
limited peaks, we can use the Single Mode Approxima- 
tion (SMA) for the dynamic correlation function. 



is') 



Jd^Q 



S{S+l) (12) 



5""(Q,cj) = S{Cl)S{huj - E{Q)), 



(10) 



for quantitative analysis of the data. Here S{Q) is the 
static structure factor, which is the Fourier transform of 
the equal-time two-spin correlation function, and quanti- 
fies the energy-integrated intensity (it is the zeroth mo- 
ment of iS""(Q,iu;)). The polarization index a may be 
omitted for an isotropic spectrum. The SMA has been 
used successfully for gapped anin chains,BO and other 
quantum many-body systemsEil where a single, coher- 
ent mode dominates the excitation spectrum. Using the 
SMA expression (Eq. in the first moment sum rule 
(Eq. ^) one can establish a relation between the measured 
intensities and the dispersion relation E(Q), 



S{Q) = - 



Sd)(l-cosQ-d) , (11) 



that involves the spin pair correlations which contribute 
to the ground state energy. 

We have carried out global fits to /m(Q, fiuj)Joi our en- 
tire data set, using Eqs. |s|-pT| as input to Eq. ^. For each 
energy scan, the non-magnetic background was modeled 
by a constant plus a Gaussian peak centered at ?ia; = 
to account for incoherent elastic scattering as needed. 
To determine the simplest parametrization of the data, 
global fits were carried out with varying numbers of terms 
included in Eq. |l^. The results of these fits with Bonds 
1-6 and Bonds 1-8 in Fig. |l| are given m Tables | and 0. 
Further neighbor bonds gave no measurable contribution 
when included in the fits. Examples of the lincshapes 
derived from the 8-bond fit arc shown as solid lines in 
Figs. ^[^. The right-hand vertical scale in these figures 
gives the normalized magnetic intensity Im (Q, fico) . Mag- 
netic neutron scattering from a local moment magnet sat- 
isfies a total moment sum rule that can be expressed as 



Clearly we do not have enough data to carry out the com- 
plete integration. However, we can ask whether the SMA, 
which accounts for the scattering data where available, 
satisfies the sum rule. Carrying out the Q-integration of 
Eq. |l| numerically, we obtain (S^) = 0.8(2). The result 
being indistinguishable from S{S + 1) = 3/4, suggests 
that the resonant mode accounted for by the SMA car- 
ries most of the spectral weight in PHCC. 

The dispersion relation determined from the 8-bond 
global fit is shown as solid lines in Figs. ^11. It is con- 
sistent with the data derived directly from the raw data, 
indicating that the variational dispersion relation Eq. y 
used in the SMA is general enough not to cause signifi- 
cant bias. As may be seen from Table the parameters 
in Eq. ^ are well determined, and the results for the 6- 
and 8-bond fits are indistinguishable. 

The first moment calculated from the model with the 
parameters determined from the global fits is shown in 
Figs. [T^p^ as dashed and solid lines for the 6-bond and 
8-bond fits, respectively. Instrumental resolution effects 
were included in this calculation. The agreement with 
the first moment determined from the Gaussian fits to 
individual scans is in general quite good, with modest 
improvement in the quality of the 8-bond fit over the 
6-bond fit. There is also good quantitative agreement 
for the terms common to both fits, particularly for the 
larger terms, indicating the robustness of the numbers 
determined through this analysis. 

Without a microscopic model that connects exchange 
constants with the dispersion relation, it is not possible 
to determine Jd and (Sq -Sd) independently. Instead one 
measures their products, which determine the contribu- 
tion of each bond to the ground state energy Eq [Eq. |j . 
In fact, because |(So • Sd)| < 3/4, each term provides a 
lower bound on the magnitude of the corresponding ex- 
change interaction. In Fig. 0(c), the thickness of the Cu- 
Cu bonds is proportional to |t/d(So-Sd)|. From this figure 
and from Table |, we see that the largest Cu-Cu,iaterac- 
tion is Bond 1, the "dimer" predicted previously .l3 How- 
ever, other bonds also show significant antiferromagnetic 
correlations, particularly Bonds 3 and 6. These bonds are 
noteworthy as they provide uniform albeit anisotropic, 
linkage of a macroscopic two dimensional spin system. 
It is also important to note that several bonds give a 
positive contribution to the first moment. These bonds, 
shown as gray lines in Fig. |l](c), arc frustrated by defini- 
tion as they raise the ground state energy. 

Summing up all terms shown in Table | we find the 
ground state energy per spin Eq — —1.5(4) meV for both 
the 6- and 8-bond fits. Using the measured value of spin 
gap A = 1 meV, we can thus determine the relation 
Eq ~ — 1.5 • A, which can be used to quantify explic- 
itly the correspondence of the spin system of PHCC to 
a particular quantum-disordered 2D spin model. It also 
helps to position PHCC on the phase diagram of the 2D 



qhafmI, characterizing how close it is to quantum crit- 
icahty. 



VI. CONCLUSIONS 

The set of measurements presented here estabhsh 
PHCC as a clean example of a two-dimensional spin- 
singlet system with a gap. The coupled-dimer picture 
that describes previously studied 2D and 3D spin gap 
materials is not appropriate for PHCC. This material lies 
in an unusual region of the phase space of the model, in a 
quantum disordered regime close to the quantum-critical 
point. Moreover, the susceptibility and heat capacity 
measurements show evidence of a quantum phase tran- 
sition from the quantum disordered to the Neel-ordered 
renormalized-classical state at a magnetic field Hd ~ 7.5 
T. Although we are not yet able to associate unambigu- 
ously the zero-field ground state with either a particular 
kind of valence bond crystal or a RVB-type spin-liquid, 
it is clear that frustration plays a key role in defining 
it. It is worth noting that evidence fop-6pin- frustration 
effects was recently observed in CuHpClE3, auather quasi- 
2D spin system, as well as in SrCu2(B03)2tll, and thus 
it may be that frustration is in fact a common feature of 
low-symmetry, gapped spin systems. 

Another issue that becomes increasingly important as 
one moves towards the QC point, or from the dimerized 
limit toward a system with more uniform coupling, is the 
appearance of a continuum in the spin fluctuation spec- 
trum. In the quantum-disordered 2D HAFM models this 
is described in terms of spinon deconfinementEj, while in 
the strongly dimerized limit it is usually adequately ac- 
counted for by multi-magnon excitations. There are some 
features in the present data indicating that a continuum 
may be detectable in PHCC [e.g. the broad, weak fea- 
ture at hu) w 3.5 meV in Fig. |^, but further experiments 
are required to clarify this issue. Finally, the low energy 
scales in PHCC make it an excellent candidate for fur- 
ther studies of quantum critical and field-dependent phe- 
nomena in a two-dimensional gapped spin system with 
frustration. 
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TABLE 1. Fractional coordinates (from Ref. [20]), bond 
lengths, and corresponding values of Jd(So • Sj) (in meV) 
for the fits of the PHCC inelastic neutron scattering data 
including 6 and 8 Cu-Cu bonds. Bond numbers correspond 
to those depicted in Fig. pl(c). 
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FIG. 1. Crystal structure of piperazinium hexachlorod- 
icuprate (PHCC), (C4Hi2N2)(Cu2Cl6). (a) View along the a 
axis showing well-separated Cu-Cl planes, (b) A single Cu-Cl 
plane viewed along the b axis, showing four unit cells. The 
dotted lines indicate possible halide-halide contacts (see text), 
(c) Cu-Cu interactions that contribute to the magnetic Hamil- 
tonian (see Table |.) The line thickness is proportional to 
|,/d(So -Sd)! for each bond. The gray lines indicate frustrated 
bonds. 



TABLE IL Fitted parameters in the dispersion relation, 
Eq. ^ for the 6- and 8-bond fits of the PHCC inelastic neutron 
scattering data. 




T(K) 

FIG. 2. Magnetic susceptibility xiT) of PHCC. Open cir- 
cles: DC suceptibility of powder. Open squares: AC sus- 
ceptibility of single crystal. Inset: solid line is a fit to a 
Curie- Weiss law. Main frame: solid line is a fit to the asymp- 
totic low-temperature susceptibility of a 2D gapped system 
as described in the text. 
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FIG. 3. AC magnetic susceptibility x{H) of single crystal 
PHCC at constant temperatures T = 0.125, 0.9, 1.2, 1.5, and 
1.8 K. Inset: Portion of H—T phase diagram of PHCC derived 
from specific heat (filled squares) and magnetic susceptibility 
(triangles) data. 
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FIG. 4. Specific heat of PHCC versus temperature at 
constant applied magnetic field, (b) shows activated behavior 
at low fields and fits used to determine field-dependence of the 
spin gap shown in (c). 
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FIG. 5. Points in reciprocal space measured using inelastic 
neutron scattering. The area of each point is proportional to 
the measured first moment of the spectrum at that Q. 
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FIG. 6. Inelastic neutron scattering data for PHCC at 
T = 1.5 K showing dispersion of magnetic excitation with 
h. Right axis shows normalized magnetic scattering intensity 
/m(Q, tj). The solid lines are a fit to a two-dimensional model 
described in the text. 
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FIG. 7. Inelastic neutron scattering data for PHCC at 
T = 1.5 K showing lack of dispersion of magnetic excitation 
with fe. Open symbols in (a): data taken at T = 50 K. Right 

axis shows normalized magnetic scattering intensity Im{^,(jj)- 
The solid lines are a fit to a two-dimensional model described 
in the text. 
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FIG. 8. Inelastic neutron scattering data for PHCC at 
T = 1.5 K showing dispersion of magnetic excitation with 
I. Right axis shows normalized magnetic scattering intensity 
Im. (Q, oj) . The solid lines are a fit to a two-dimensional model 
described in the text. 
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FIG. 9. Energy of magnetic excitations in PHCC showing 
dispersion with h at constant k and I. Data points are deter- 
mined from Gaussian fits to constant-Q scans. The solid Unes 
are a fit to a two-dimensional model described in the text. 
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FIG. 10. Energy of magnetic excitations in PHCC show- 
ing lack of dispersion with k at constant h and I. Data points 
are determined from Gaussian fits to constant-Q scans. The 
solid lines are fit to a two-dimensional model described in the 
text. 




FIG. 12. Variation of first moment h{uj)Q of PHCC with 
FIG. 11. Energy of magnetic excitations in PHCC show- h at constant k and I. Data points are determined by fitting 
ing dispersion with I at constant h and k. Data points arc c^ch original constant-Q scan to a Gaussian peak and calcu- 
determined from Gaussian fits to constant-Q scans. The soUd lating the first moment of the intensity from that fit. The 
Unes are fit to a two-dimensional model described in the text. dashed (solid) Unes are determined by a global fit including 



6(8) Cu-Cu interactions, as described in the text. 
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FIG. 13. Variation of first moment Tl{u))ci of PHCC with 
k at constant /i and I. Data points are determined by fitting 
each original constant-Q scan to a Gaussian peak and calcu- 
lating the first moment of the intensity from that fit. The 
dashed (soUd) lines are determined by a global fit including 
6(8) Cu-Cu interactions, as described in the text. 
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FIG. 14. Variation of first moment ?i(w)q of PHCC with 
I at constant h and k. Data points are determined by fitting 
each original constant-Q scan to a Gaussian peak and calcu- 
lating the first moment of the intensity from that fit. The 
dashed (solid) lines are determined by a global fit including 
6(8) Cu-Cu interactions, as described in the text. 



